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In this paper we investigate the role of the electronic correlation on the hole doping dependence 
of electron-phonon and superconducting properties of cuprates. We introduce a simple analytical 
expression for the one-particle Green's function in the presence of electronic correlation and we eval- 
uate the reduction of the screening properties as the electronic correlation increases by approaching 
half-filling. The poor screening properties play an important role within the context of the nona- 
diabatic theory of superconductivity. We show that a consistent inclusion of the reduced screening 
properties in the nonadiabatic theory can account in a natural way for the Tc-S phase diagram of 
cuprates. Experimental evidences are also discussed. 
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I. INTRODUCTION 

The role of the electron-phonon (el-ph) interaction in 
the high- Tc superconducting cuprates has been a mat- 
ter of debate for a long time. In early times the report 
of a negligible isotope effect on Tc at optimal doping, 
the almost linear behaviour of the resistivity on tem- 
perature, also at optimal doping, and other exotic fea- 
tures of the copper oxides led to the common belief that 
electron-phonon coupling was a marginal ingredient to 
understand the phenomenology of these materials. How- 
ever, over the years, there has been a revamping evidence 
of an important role of the phonons. The most remark- 
able ones are, for instance, the discovery of an isotope 
effect on Tc larger than the BCS value {ar^ > 0.5) in 
the undcrdoped regime,J^ the report of a sizable iso- 
tope shift on the effective electronic mass m*— and on 
the onset of the pseudogapr^ the observation of phonon 
renormalization^ and phonon anomalies at T < Tci 
More recently, ARPES measurements pointed out a kink 
in the electron dispersion the origin of which is probably 
phononici^ Clearly, if phonons are relevant for supercon- 
ductivity in these materials, this cannot be described in 
a BCS-like framework, but some non-conventional ap- 
proach including strong electronic correlation is neces- 
sary. The study of the interplay between electron-phonon 
interaction and the electronic correlation is a challenging 
task which has attracted much scientific work along dif- 
ferent lines. 

An interesting issue concerns the momentum modu- 
lation of the electron-phonon coupling induced by the 
electronic correlation. In Ref. 0, using a variety of theo- 
retical and experimental findings, it is shown that in cor- 
related systems small-q scattering in the electron-phonon 
interaction is strongly favored. A strong enhancement of 
the forward scattering at q ~ in correlated systems 
close to the metal-insulator transition, accompanied by a 
suppression of scatte ring at large q, was reported for ex- 
ample in Refs. Ildlll bv using expansion techniques. 



A recent numerical work based on Quantum Monte Carlo 
technique confirms this picture 

Different but somehow complementary argumentations 
based on poor screening effects in correlated systems have 
been also discussed in literature. The basic idea is that, 
as a metal loses its coherence as function of the corre- 
lation degree approaching a metal-insulator transition, 
the screening properties of the bare long-range electron- 
phonon interaction become less effective resulting in a net 
predominance of small q scatteringii2iiiiiS4i& A similar 
physical argument applies, for example, to doped semi- 
conductors which are commonly described in terms of 
the Frohlich Hamiltonian, with electron-phonon matrix 
elements l^qp cx l/|qp. 

The momentum structure of the electron-phonon scat- 
tering induced by the electronic correlation has been 
shown to play a crucial role in the context of nonadia- 
batic superconductivity^iiSiiS In narrow band systems, 
such as cuprates and fuUerenes, the Fermi energy Ep is 
so small to be comparable with the phonon frequencies 
tjph, and the adiabatic assumption (cjph ^ Ep) breaks 
down. In this context Migdal's theorem.^° does not apply 
and one needs to take into account nonadiabatic effects 
not included in the Migdal-Eliashberg (ME) theory of 
superconductivity. Detailed studies have shown that the 
nonadiabatic contributions, which are well represented 
by the vertex function, present a complex momentum- 
frequency structure, in which small q-scattering leads to 
an enhancement of the effective superconducting pairing, 
while large-q scattering leads to a reduction of itiiiiSiiS 
The strong q-modulation of the electron-phonon interac- 
tion due to the electronic correlation is thus expected to 
give rise to a net enhancement of the superconducting 
pairing. 

The purpose of the present paper is twofold. On one 
hand we wish to quantify the microscopic dependence 
of the screening properties of a correlated system on 
relevant quantities as the electron density of the Hub- 
bard repulsion; in addition we apply the derived screened 
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electron-phonon interaction to evaluate the role of the 
electronic correlation in the context of the nonadiabatic 
superconductivity and to derive a qualitative supercon- 
ducting phase diagram. To this aim we introduce a 
model for the electronic Green's function of the sys- 
tem, based on the decomposition of the total spectral 
function in a coherent, itinerant part, and an incoher- 
ent localized background corresponding to the Hubbard 
subbands. The relative balance between the two parts 
varies as a function of doping and electronic correlation. 
This will have important consequences on the electronic 
screening and hence on the q-modulation of the effective 
electron-phonon scattering, as well as on the supercon- 
ducting properties. We shall show that: 

• the coherent excitations dominate the screening 
properties as well as the superconducting ones. 

• the loss of coherent spectral weight approaching 
half-filling is thus responsible for the reduction of 
the screening properties and for the increase of the 
forward scattering in the electron-phonon interac- 
tion. 

• in the strongly correlated regime the selection of 
forward scattering gives rise to an enhancement 
of the effective electron-phonon interaction within 
the context of the nonadiabatic superconductivity. 
These effects however compete with the reduction 
of the quasi-particle spectral weight which is detri- 
mental for superconductivity. 

• the resulting phase diagram shares many similari- 
ties with the one of the cuprates. In particular it 
shows an overdoped region, where superconductiv- 
ity is suppressed by negative nonadiabatic effects, 
an underdoped region, in which superconductiv- 
ity is destroyed by the loss of coherent spectral 
weight, and an intermediate region in which the 
predominance of small-q scattering leads to an en- 
hancement of the nonadiabatic el-ph pairing which 
overcomes the reduction of the coherent spectral 
weight. 

We hereby wish to point out that a complete descrip- 
tion of the rich fenomenology of the cuprates is well be- 
yond the aim of the present paper. In particular, we 
shall not discuss, for reason of simplicity, the symmetry 
of the gap, which of course is of fundamental importance 
if one wishes to give a quantitative description of these 
systems. We would like just to remark on this point that 
a d-wave symmetry of the superconducting order param- 
eter was shown by many authors to naturally arise in the 
context of a phonon pairing with a significant predomi- 
nance of forward scattering2ii2SiSiS4. The competition 
between s— and d— wave symmetry in a nonadiabatic 
electron-phonon system was also studied in Ref. Tak- 
ing into account explicitely the d-wave symmetry of the 
gap would not change in a qualitative way the results of 
the present work. 



This paper is organized as follows: in section ^ we 
introduce our model Green's function; in section ITTTl we 
derive an effective form for the electron-phonon interac- 
tion. In the last section we write and solve the generalized 
Migdal-Eliashberg equations, in the adiabatic and non- 
adiabatic limit, and discuss in detail the competition of 
the different factors which determine the superconduct- 
ing critical temperature of our system. 

II. A MODEL FOR CORRELATED ELECTRON 
SYSTEMS 

As briefly discussed in the introduction, one of the 
main aims of the present paper is to investigate how the 
screening properties are affected by the presence of strong 
electronic correlation, and to parametrize these effects in 
terms of microscopical quantities. In particular we have 
in mind a Hubbard-like system where itinerant electrons, 
with band dispersion and bandwidth i?, interact each 
other through an onsite Coulomb repulsion U. As we are 
going to see, a crucial role is played in this context by 
the transfer of spectral weight as a function of the corre- 
lation degree from low energy coherent states to the high 
energy (Hubbard-like) incoherent ones. 

In this section we present a simple, minimal model for 
the electron spectral function which takes into account 
these main effects and which can thus represent a proper 
starting point to evaluate screening effects in correlated 
systems. 

All the possible information about the single-particle 
properties of the system is contained in the one-electron 
Green's function G(k, lj). Without loss of generality we 
assume that the Green's function G can be split in a 
coherent and an incoherent contribution^^ 

G(k, iu) = Gcoh(k, uj) + Gi„c(k, u;), (1) 

where the coherent part Gcoh describes the itinerant, 
quasi-particle like properties of the electron wavefunc- 
tion, while the incoherent part Ginc accounts for the in- 
coherent high energy excitations. Due to its localized 
nature Ginc(k,w) is only weakly dependent on the mo- 
mentum quantum number, so that the dependence on k 
can be reasonably neglected. 

An important quantity which parametrizes the relative 
balance between coherent and incoherent contributions 
is the quasi-particle spectral weight Z, which is simply 
given by: 

J diuhm [Gcoh(k, UJ + iS)] = Z, (2) 

whereas the incoherent part obeys the sum rule: 

J dcj^Im [Gi„c(k, UJ + iS)] = 1- Z. (3) 

The quasi-particle spectral weight Z can vary between 
and 1, the two limits corresponding to the insulating 
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and metallic limit respectively. It depends on the internal 
parameters U and S, where S is the hole doping (S = 1—n) 
and n the total number of electrons (n = 1 half-filled 
case). 

Several techniques have been developed to investigate 
the Hubbard modeliSl Different starting points are em- 
ployed according to whether main emphasis has to be 
paid on the coherent (itinerant) or on the incoherent (lo- 
calized) features. For instance the so-called Hubbard I 
approximatioufS which is exact in the atomic limit, is 
mainly aimed at a schematic representation of the local- 
ized states, described by an upper and a lower Hubbard 
band spaced by an energy gap of width U. On the other 
hand the Gutzwiller technique^S and the mean field slave 
bosons solution^" offer an useful tool to deal with the co- 
herent spectral weight of the electron Green's function: 
in this case the quasi-particle spectral properties in the 
presence of strong correlation are described in terms of an 
effective band of non-interacting fermions with spectral 
weight Z and bandwidth ZE. 

In this paper we introduce a new phenomenological 
model to take into account in the simplest way and at 
the same level the coherent and incoherent parts of the 
Green's function. We approximate the exact (unknown) 
coherent and incoherent parts of G(k, oj) in Eq. ^ re- 
spectively with the Gutzwiller— and Hubbard solu- 
tions, namely: 



G'coh(k, w) 



(1-71/2) 



(4) 



w - (1 - n/2)ek + /i - C//2 
n/2 



w - (n/2)ek + /i + 



(5) 



where fi is the chemical potential, Ng is the total num- 
ber of sites and Z is the quasi-particle weight obtained in 
the Gutzwiller approximation in the paramagnetic state 
at finite U and generic filling (Appendix Due to 

the localized nature of the incoherent part we have re- 
placed the Ginc (k, to) given by the Hubbard I approxima- 
tion with its momentum average. Numerical calculations 
based on Dynamical Mean-Field Theory (DMFT) con- 
firm our qualitative picture of a spectral weight transfer 
from a central coherent peak to a incoherent Hubbard- 
like background with increasing U 

The behaviour of Z as function of the particles den- 
sity n and of the Hubbard energy U is shown in Fig. 
^ The critical Hubbard energy Uc, which determines 
the Brinkman-Rice transition at n = 1 is related to the 
kinetic energy i?kin, which depends on the bare electron 
dispersion shape, through the relation Uc = 8|i?kin| ^ In 
the following we employ a bare constant density of states 
(DOS) with 7V(ek) ^ Nq ^ 1/E for et G [-E/2,E/2]. 
In this case, we have Uc = 2E. The chemical potential /i 
is determined by the total number of particles. In Fig. [21 
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FIG. 1: Quasi-particle spectral weight as determined by 
the Gutzwiller solution at finite U and n. Left panel: Z 
as function of n for (from top to the bottom) U/Uc = 
0.4,0.8, 1.2, 1.6,2.0. Right panel: Z as function of U/Uc for 
(from top to the bottom) n = 0.6, 0.8, 0.9, 1.0. 



we show typical density of states N{uj) for the correlated 
system described by our model [Eqs. Q-ljSJ]. 

We would like to stress that the phenomenological 
model described by Eqs. Q-© is not meant at all to 
be exhaustive of the complex physics of a strongly cor- 
related system. In fact, retardation effects are neglected, 
since we are assuming the separation into two species of 
electrons to be independent of frequency. More sofisti- 
cated methods of solution, including DMFT, permit to 
treat the self-energy of a strongly correlated system in a 
more careful way, retaining the correct frequency depen- 
dence of the self-energy. 

Our model has the advantage of being extremely sim- 
ple and easy to handle, and it allowed us to obtain ex- 
plicit expressions for all the relevant quantities of the 
coupled electron-phonon system; in particular, we focus 
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FIG. 2: Density of states N{uj) = (I/tt) ^^^'(k, (.j -I- i0+) 
resulting from our model, for U — 2Uc and different values of 
doping. At half-filling the system is an insulator, and its den- 
sity of states is represented by two Hubbard-like features at 
distance U from each other; moving away from half-filling 
a coherent peak starts forming, with increasing weight Z. 
Dashed regions represent filled states up to the chemical po- 
tential /i (dotted line). 
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on the spectral weight transfer from the coherent to the 
incoherent part of the Green's function when increasing 
the degree of electronic correlation. As we are going to 
see, this feature will have important consequences on the 
electronic screening and on the momentum dependence 
of the electron-phonon coupling. 



III. SCREENING AND ELECTRON-PHONON 
INTERACTION 

A. Correlation effects on Thomas-Fermi screening 

The momentum dependence of the electron-phonon in- 
teraction usually plays a marginal role in determining the 
electron-phonon properties of common metals. The basi- 
lar reason for this is that the bare long-range electron- 
phonon interaction is effectively screened by the long- 
range Coulomb repulsion leading to a weak momentum 
dependence. 

The conventional Migdal-Eliashberg theory, which de- 
scribes electron-phonon effects both of the normal and 
superconducting states, is formally derived starting from 
an effective electron-phonon Hamiltonian, in which the 
Coulomb electron-electron repulsion does not appear, 
apart from a weak residual electron-electron contribution 
in the Cooper channel, C/k,q, which gives rise the to the 
Morel- Anderson "pseudopotential" term /x = N{0)U 
The physical quantities appearing in this effective Hamil- 
tonian are thus considered to have been already renor- 
malized by the long-range Coulomb interaction. In par- 
ticular the electron-phonon matrix elements 5k.k+q and 
the residual electron-electron repulsion are usually con- 
sidered to have a negligible momentum dependence, so 
that the Eliashberg equations depend only on the fre- 
quency variables. 

This drastic assumption works quite well in many con- 
ventional low temperature superconductors with large 
carrier density since, in this case, the long-range q- 
dependence of the bare electron-phonon and electron- 
electron interaction [^(q, uj) oc l/|qp)] is removed by the 
large metallic screening. This well-known effect is usu- 
ally expressed in terms of the (static) dielectric function 
e(q), which in the RPA approximation reads: 



e(q) = l 



(6) 



where fcxF is the Thomas- Fermi screening momentum de- 
fined as 



fc|p = - lim 47re^n(q, w = 0), 

q^O 



(7) 



and n(q, uj): 

n(q, u;)^^Y. [ ^'^'^(l^ + q> ^ + ^')G(k, Lu'). (8) 



The effective long-range interaction results thus screened 
by conduction charge to give the Thomas- Fermi expres- 



Foff(q,t^) = 



cx 



<q) 
1 

Iqp + fcTp' 



(9) 



In free electron systems the Thomas-Fermi vector is 
directly related to the bare density of states via the sim- 
ple relation limq^on(q,tj = 0) = -2iV(0), where iV(0) 
is the density of states per spin at the Fermi level, so 
that fc|p = 87re27V(0). In common metals, since fcxF is 
typically larger than the Brillouin zone size (fcez), the 
effective (electron-electron, electron-phonon) interaction 
Vcff{(l, to) can be considered in first approximation almost 
independent of the exchanged momentum q. 

Things are expected to be very different in correlated, 
narrow band systems. As we have mentioned before, 
strongly correlated electrons, due to their reduced mo- 
bility, are much less effective in screening external per- 
turbations, especially at small wavelengths. For instance, 
the reduction of the screening properties approaching a 
metal-insulator transition in disorder alloys as well as in 
cuprates has been experimentally signaled in Refs. I40ll4ll 

In this section we employ the simple model above in- 
troduced for the description of the Green's function to 
quantify the reduced screening properties of correlated 
systems and their dependence on microscopic parame- 
ters, such as the hole doping i5 or the Hubbard repulsion 
U. In order to do this, we compute the Thomas-Fermi 
vector fcxF , defined in Eq. Q , using the model described 
by Eqs. Q-© to evaluate the RPA response function 
n(q, uj) according to Eq. JSJ. While higher order (vertex) 
diagrams are not taken into account in this framework, 
we shall show that this simple model is already sufficient 
to describe the reduction of screening properties due to 
transfer of spectral weight from the coherent to incoher- 
ent states. 

Using Eqs. Q-© the response function H can be 
written as a sum of three different contributions: 



H = n, 



Hi- 



(10) 



where the first one describes scattering processes which 
involve only coherent states; the second term describes 
scattering between the coherent peak and the Hubbard 
lower/upper (incoherent) bands; the last one describes 
processes which involve only localized incoherent states 
in both the Green's functions of Eq. © . In general we 
expect that the total screening will be dominated by the 
first contribution Hc_c since the itinerant coherent states 
are much more effective, because of their mobility, in 
screening external perturbations than the localized ones. 

In Fig. Owe plot the RPA response function in units of 
the bare DOS: — limq^on(q, w = 0)/27Vo7 as function of 
the electron filling. Since for n — s- the screening proper- 
ties are determined only by the coherent part regardless 
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FIG. 3: Effective Thomas-Fermi screening fc^p (solid line) as 
a function of the electron density n for a correlated system 
described by our model {U = SUc)- The different contribu- 
tions to the total screening are also shown: coherent-coherent 
particle-hole processes (dashed line) and coherent-incoherent 
+ incoherent-incoherent contribution (dotted line). Apart 
from half-filling, where the coherent contribution vanishes and 
the screening is determined by the only residual incoherent 
polarization, the static screening properties of the system are 
dominated from the coherent quasi-particle excitations. 



any correlation effects, this is also equivalent to plot the 
Thomas-Fermi momentum k^p as function of the elec- 
tron density n: k^p(n)/k^p{n = 0). As shown in figure 
the net value of the Thomas- Fermi momentum is mainly 
determined by the coherent-coherent excitations. Simple 
scaling considerations show that the coherent-coherent 
contribution to the response function is just equal to 
times the Thomas-Fermi momentum of an uncorrelated 
system. The explicit expressions of the other two terms 
are a bit more complicated and they are reported in ap- 
pendix Fig. 121 shows a drastic reduction of the screen- 
ing properties of the system as the metal-insulator tran- 
sition is approached at half-filhng {U > Uc}- In this case 
the spectral weight of the coherent part is zero, and the 
only residual small contribution to the screening is due 
to incoherent excitations which vanishes for U oo. 



B. Poor screening and momentum dependence of 
the electron-phonon interaction 

From Fig. |2| it is clear that the assumption of a 
Thomas-Fermi momentum much larger than the ex- 
changed phonon momenta q breaks down as electronic 
correlation effects get more and more relevant, namely 
approaching half-filling. In this situation the effective 
electron-phonon interaction can be no longer considered 
weakly dependent on q in the long-range limit q — > 0. 
On a microscopical ground the screening of long-range 



Coulomb interaction renormalizes both the bare electron- 
phonon matrix element ^^^^ and the phonon frequen- 
cies riq. The el-ph matrix element can be usefully writ- 
ten as ~ c(riq)/|q|, where c(riq) is a well behaved 
function of q in the limit limq^o and it mainly depends 
on the phonon frequency Qq. If both the screening ef- 
fects on g{J and flq are properly taken into account^ 
one can get an expression for the effective total electron- 
phonon interaction: 

= (11) 

where both the phonon propagator Dq{Lo) and the cou- 
pling function c(a;q) are written in terms of the screened 
phonon frequency Wq. Eq. shows that the long-range 
behaviour of the total el-ph interaction oc l/|qp, when 
written as function of the screened phonon frequency, is 
correct by the dielectric function e(q). 

For an optical mode, Wq is only weakly dependent on q 
and the leading dependence on q of Eq. comes from 
the term cx l/[e(q)|qp]. These screening effects can be 
conveniently dealt with by introducing the screened el-ph 
matrix element ^q: 

2 li'ql ^ /ION 

The el-ph scattering is thus roughly described (we remind 
these expression were derived in the limit q — ^ 0) by a 
lorentzian function in the space |q|. It is also useful to 
introduce the dimensionless variables Q — |q|/2fcF and 
Qc = ^Tp/Sfep, so that: 

The parameter Qc represents a cut-off for the exchanged 
phonon momenta: the electron-phonon scattering will be 
operative for Q ^ Qc, and negligible for Q ^ Qc- 

The momentum structure resulting in Eq. plays 
a crucial role in the Cooper pairing in the coherent- 
coherent channel where the momentum is a good quan- 
tum number. For these contributions the total strength 
of the electron-phonon coupling is linked with the mo- 
mentum average of Eq. over the Fermi surface. 
For a isotropic system, using polar coordinates J dfl — 
J^^ d(f) dcosO and reminding that Q = sin{6/2), we 
obtain: 

Jdcl^JoQdQ-^^j^ 

In common metals Qc ~ 0.5 — 1 so that In [(1 -I- Qc)/Qc] 
is of the order of 1. On the other hand, in poorly screened 
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systems Qc ^ 1 and the resulting el-ph coupling is sensi- 
bly enhanced. In the following we shall consider — 0.7 
as representative case of uncorrelated usu al ineta ls. 

For practical purposes, following Refs. 'iS'lO", we ap- 
proximate the lorentzian behaviour of Eq. H13|l with a 
Heaviside 9 function: 

\9{Q)? ^ g^^mc-Q). (15) 

In order to preserve in this mapping the total strength 
of the el-ph coupling, the prefactor a has to be deter- 
mined by requiring the resulting el-ph coupling strength, 
namely the average of over the Fermi surface, to be 
equal for Eqs. H13() and (|15() . With this condition we 
find: 

= In (^^) ^(Q^ - Q)' (16) 

As a final remark of this section we note that the mo- 
mentum dependence of |5(Q)P is not expected on the 
other hand to be effective in the incoherent-coherent 
and incoherent-incoherent contributions to the electron- 
phonon interaction, where the exchanged momentum q 
is no more a good quantum number. In this case the 
effective incoherent electron-phonon coupling is roughly 
given by its momentum average on the Brillouin zone, 
which we shall set in the following to be equal to . 

IV. GENERALIZED MIGDAL-ELIASHBERG 
EQUATIONS 

In the previous sections we have introduced a simple 
model for an electron-phonon system in the presence of 
electronic correlation. In particular we have reduced, in 
an approximate way, the complex problem of the inter- 
play between electron-phonon and electron-electron in- 
teractions to a purely electron-phonon system described 
by an effective one-particle Green's function [Eqs. iQ), 
Q , ©J and an effective electron-phonon matrix element 
g{Q) [Eq. H16|l ]. After this mapping, the Baym-Kadanoff 
theory— ^ assures that the functional form of the supercon- 
ducting equations will be the same of a purely electron- 
phonon system: 

$ = $el-ph[<7,G,$], (17) 

Z = Zei-ph[5,G,$], (18) 

where $ is the superconducting order parameter; the 
Green's function G and the matrix element g are de- 
fined by Eqs. ((SJ, H16|l . as mentioned above. In 
order to obtain an explicit expression for Eqs. H17|l - H18|) 
we should specify in which framework we are going to 
treat the electron phonon interaction. In particular, we 
observe that the conventional ME theory, in particular, is 
based on the assumption that the phonon frequencies are 
much smaller than the electronic Fermi energy, ^ 
( adiabatic limit). This theory works quite well in the 



conventional low temperature superconductors, where no 
electronic correlation is present and E-p is of the order of 
5—10 eV. On the other hand, the strong band renormal- 
ization in correlated systems described in Sec. ^ ques- 
tions the adiabatic assumption, especially as, approach- 
ing half-filling, the renormalized bandwidth ^ ZE can be 
comparable with Wph. In these systems a more suitable 
description can be obtained in the framework of the non 
adiabatic theory of superconductivity>i2ii2iiS Eqs. 117() - 
(|18|l can be rewritten as: 

Zn = 1 + — Vrz([G];w„,c^„)r;„[G], (19) 

— 

m 

<&„ = T,^r$([G];c.„,c.„)^^,7^^[G], (20) 

m 

where the electron-phonon kernels rz([G];w„,a;,„) and 
r$([G]; a;„, cjm) contain the nonadiabatic vertex (P) and 
cross (G) contributions to the self-energy and to the 
Cooper pairing channels: 

rz([G];LJ„,cj„i) = A„_,„ [1 + AP([G];t^„, Qc)] , 
r*([G]; [1 + 2AP([G]; w„,a;,„, Qc)] 

+A^G([G]; w„,Ci;„i, Qc) - /i. 

Here Xn-m is linked with the electron-phonon spec- 
tral function a^F{ijj) through the relation A„_m = 
2 /dwa^F(w)w/[w2 -I- [bJn - ujm)% A = A„_m=o and 

is the short-range residual electron-electron repul- 
sion. The breakdown of the adiabatic hypothesis de- 
termines the need for the explicit inclusion of the ver- 
tex (P) and cross (G) functions in Eqs. ((T^ - lf^ and 
it affects the expression of /^^[G] = X]k^(^''^) ^^"^ 
vt[G] = EkG(k,w)G(-k,-u;) through finite band- 
width effects. The momentum dependence of the super- 
conducting equations has been averaged on the Fermi 
surface and it gives rise to the strong dependence on Qc 
in the vertex and cross terms. In Eq. 119|) - (|2U|I we have 
moreover implicitly expressed the functional dependence 
of the electron-phonon kernels F^, Fa as well as of the 
quantities P, G and 77, on the Green's function G which 
we remind is modeled as in Eqs. QJ, (jSJ. 

Before solving Eqs. (|19|I - H20|I in the whole range of dop- 
ing, we would like to discuss the different role of the co- 
herent (itinerant) states and the incoherent (localized) 
states, described respectively by Eqs. © on the su- 
perconducting properties. As we have seen in Sec. IIIII 
the electronic screening is mainly dominated by the co- 
herent term of the electronic Green's function, which de- 
scribes mobile electrons for which k is a good quantum 
number. 

Similar considerations can be made also for supercon- 
ductivity: we expect, in fact, that the coherent electrons, 
which have a high mobility, will give a more relevant 
contribution to the superconducting critical temperature. 
To check the validity of this hypothesis, we have solved 
Eqs. lfTO|l - lf^ in the ME limit (i.e., neglecting vertex 
corrections) , once using an integral kernel containing the 
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0.25 




u/u, 

c 

FIG. 4: Comparison of the critical temperature Tc as a func- 
tion of U for the half-filling case, using the full integral kernel 
in Eqs. 11911 - 120(1 (empty circles) and using only its coherent 
part (solid line). 



full Green's function (Eq. and once an integral ker- 
nel with only the coherent part of the Green's function 
[Eq. Q], as a function of the Hubbard repulsion U . In 
Fig. we show as empty circles the results obtained 
with the full kernel, and with solid line the critical tem- 
perature obtained using only the coherent part. The two 
sets of data are hardly distinguishable, pointing out that 
the increase of Tc due to the coherent-incoherent and 
incoherent-incoherent couplings is negligible. 

After this observation, in the following the functional 
dependence on the total Green's function G in Eqs. H19() - 
(I2()|l can be in good approximation replaced by the only 
coherent part, explicitly: T z[G coh\, r$[G'coh], P[G co\y], 
C[Gcoh], 'nm\G coh\- As wc show in Appendix 151 when the 
reduced spectral weight and bandwidth are taken into 
account, this corresponds to a proper rescaling of the 
analytical expressions for these q uantiti es evaluated in 
the absence of correlation in Refs. 

It is interesting to compare our model with the two- 
band superconductivity^^, which has recently driven a 
considerable attention due to MgB2« In that case, 
the opening of inter-band scattering channels leads to 
an enhancement of the critical temperature. For some 
respects, our model could also be seen as an effective 
two-band system, made up of a very narrow band of mo- 
bile electrons and another band of localized electrons, 
coupled to each other. However, we note that, since the 
spectral weight of each single band is not conserved, the 
onset of the high-energy bands of localized electrons is 
accompanied by a decrease of the quasi-particle spectral 
weight, resulting in an effective reduction of the Cooper 
pairing. 



A. Doping effects and phase diagram of the 
nonadiabatic superconductivity 

Eqs. lfTn|) - (pn|l represent our tool to investigate the 
loss of the superconducting properties due to the elec- 
tronic correlation approaching half-filling. We can in fact 
evaluate all the relevant quantities, such as the electron- 
phonon interaction kernels F^, Fa, the electron Green's 
function G, the vertex and cross functions P, C, and the 
momentum cut-off Qc as a function of the microscopic pa- 
rameters as the hole doping 6 and the Hubbard repulsion 
U. We shall show that the phase diagram as a function 
of the doping is governed by two competing effects: one 
driven by the reduction of the coherent spectral weight 
approaching half-filling, which is detrimental for super- 
conductivity, and the other by the complex behaviour 
of the non-adiabatic terms, which increase the effective 
pairing as 5 — s- and decrease it as (5 —> 1. 

Since we are mainly interested in the region (5 — > of 
the phase diagram, we disregard for simplicity the an- 
alytical dependence of the non-adiabatic terms on the 
chemical potential. The behaviour of the "bare" P and 
C as a function of doping is in fact determined by the 
density of electrons {n = 1 — S); this dependence is much 
weaker than the dependence of Z and Qc close to half- 
filling (see Figs. El III). 

Before solving Eqs. (|19|I - H2U|) numerically to obtain 
the critical temperature Tc as a function of doping, we 
wish to discuss the phase diagram of our model in terms 
of simple intuitive physical arguments, based on an ef- 
fective electron-phonon coupling. Let us consider for the 
moment the electron-phonon interaction alone, without 
any residual Coulomb repulsion, namely ^ — 0. Eq. H20() 
can be rewritten in a simplified way as: 

$„ ~ Te ^ ZA [1 + 2ZXP + ZXC] K^-,n^m, (21) 
m 

where we have simplified, according Appendix ^ the 
main dependences on Z = Z{U, S) of each quantity. 
In this way, we can roughly see the total electron- 
phonon coupling as the product of two terms: an effective 
electron-phonon coupling of ME theory renormalized by 
the electronic correlation, A^^, and the enhancement due 
to nonadiabatic vertex and cross (VC) diagrams 7^'-': 

A<=« = A^E7^'', 
A^^ = ZA, 

7VC ^ l + 2ZXP{Qc) + ZXCiQc). 

The schematic behaviour of these quantities as a function 
of the hole doping S is shown in the upper panel of Fig. |S1 
The physics behind the (5-dependence of A^^ can be eas- 
ily related to the loss of spectral weight approaching the 
metal-insulator transition for 5 ^ 0. This effect, which 
is present also in 7^*^^, is however in that case compet- 
ing with the enhancement of the effective coupling due 
to P(Qc) and C{Qc) which will be maximum and posi- 
tive close to half-filling (where Qc 0) and negative at 
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FIG. 5: Graphical sketch of the difFerent contributions to the 
effective superconducting couphng. Top paneh the couphng 
function is mainly determined by the coherent spectral 
weight, and it exhibits a monotonous growing behaviour as a 
function of doping. The vertex factor 7^^^ tends to enhance 
the effective coupling at low doping and to depress it at high 
doping. Middle panel: the total effective electron-phonon cou- 
pling A"''^ = x'^^y^'^ has a maximum at some finite value 
of S; when the effective Morel- Anderson pseudopotential is 
subtracted, superconductivity is suppressed at high doping. 
Lower panel: resulting phase diagram for superconductivity: 
superconductivity is only possible in a finite region of phase 
space (gray region), where A^*^ 



fi'^'^ is positive. 



high dopings. The interplay between these two effects 
will give rise to a maximum of 7^^, and hence of A°*, 
somewhere in the small doping region where the compe- 
tition between the spectral weight loss and the positive 
nonadiabatic effects is stronger (see upper and middle 
panels in Fig. 

We can now also consider the effect of the residual 
Morel- Anderson-like repulsion /x; first of all, we observe 
that the reduction of spectral weight will lead to an ef- 
fective repulsion /i^*^ ~ Z^. Superconductivity will be 
possible only when the net electron-phonon attraction 
overcomes the repulsion term: — > (see lower 
panel of Fig. Isj. The resulting total coupling is ex- 
pected to exhibit a "bell" shape which is mostly due to 
the (5-dependence of the nonadiabatic factor 7^*-^. It is 
interesting to note two things. First, in the extreme case 




FIG. 6: Maximum eigenvalue as a function of doping, 

evaluated at T = O.Oloio for different values of fi, with A = 
1, uJo = 0.8-E/2, and U = 8Uc. Empty symbols (dashed 
hues) represent ME theory, filled symbols (solid lines) the 
nonadiabatic theory described by Eqs. I191 - I2m . 



^ME < ^off ^ where no superconductivity would be pre- 
dicted in the whole S range by the conventional ME the- 
ory, we could expect finite Tc in a small 6 region, due to 
purely nonadiabatic effects A°*^ — A'^^7^*-^ > . Sec- 
ondly, it is clear that within the ME framework a net 
attractive interaction in the Cooper channel at a certain 



doping 5, which corresponds to A*^^ > , would imply 
a superconducting order also at larger S since the two 
quantities X^^,ii''^ scale in the same way cx Z; on the 
other hand, in the nonadiabatic theory superconductiv- 
ity, Tc is expected to be limited to some maximum value 
of doping, due to the negative contribution of the nona- 
diabatic diagrams P and C at large 6 (large Qcs). 

We can now quantify the simple arguments discussed 
so far. A quantitative estimate of the strength of the 
superconducting pairing is given by the highest eigen- 
value v"'^'^^ of the superconducting integral kernel in Eq. 
(|20|l . computed at low T; at a given temperature T and 
doping S superconductivity occurs if v"^^^ > 1 and the 
superconducting pairing (and Tc) is stronger as v"^'^^ is 
larger. 

In Fig. El we compare the behaviour of v''^^^ as a bmc- 
tion of 6, obtained at T = O.OIwq using an Einstein 
spectrum for different values of fj, in ME (open sym- 
bols, dashed lines) and in the nonadiabatic theory (full 
symbols, solid lines). The Hubbard repulsion was set at 
U — 8Uc and the phonon frequency at ujq — 0.8E/2, 
where E/2 is the bare half-bandwidth (unrenormalized 
by correlation effects) . The corresponding phase diagram 
Tc vs. 5 is reported in Fig. In agreement with our pre- 
vious discussion in the ME framework y'^'^^ and Tr de- 
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FIG. 7: Phase diagram of Tc vs. 5 for the adiabatic (empty 
symbols, dashed hnes) and nonadiabatic (fiUed symbols, solid 
lines) theory. Details as in Fig. |S| 



crease monotonously as the hole doping is reduced. On 
the other hand, the corresponding results in the nonadia- 
batic theory display a more complex behaviour, showing 
that the effective nonadiabatic pairing is larger than the 
ME one at low doping and smaller at high doping. 

As we have already discussed, the bell-shape of the 
highest eigenvalue v"'^'^^ and of the critical temperature 
Tc can be attributed to the dependence of the magnitude 
and sign of the nonadiabatic terms on Qc, which, in turn, 
strongly depends on doping. For high doping the nonadi- 
abatic contributions are negative and decrease ■y™'*'^ and 
Tc with respect to their ME values. Decreasing S the 
nonadiabatic terms turn from negative into positive and 
^max g^j-^^ increase up to a maximum value. As the hole 
doping is further decreased ((5^0), the loss of spectral 
weight becomes the dominant effect and it finally leads to 
the complete suppression of superconductivity. The in- 
clusion of the residual Coulomb repulsion in the Cooper 
channel, fi, leads to an overall reduction of the supercon- 
ducting pairing. The effect is more pronounced in the 
nonadiabatic theory than in ME, since in this case a very 
small value of is enough to suppress superconductivity 
in a large region of phase space at high dopings. 



V. DISCUSSION AND CONCLUSIONS 

The main aim of the present paper is the description, 
on microscopical grounds, of the hole doping dependence 
of the electron-phonon and superconducting properties of 
a strongly correlated system within the context of a nona- 
diabatic electron-phonon theory. The need for a nonadi- 
abatic treatment of the electron-phonon interaction in 
correlated systems comes from the fact that, as a metal- 



insulator transition is approached, the electronic band- 
width is strongly reduced, and the adiabatic assumption 
i^ph/Ep on which Migdal's theorem is based breaks down. 

Past studies have shown that the inclusion of nonadi- 
abatic effects can lead to a strong enhancement or de- 
pression of Tc depending on the value of the exchanged 
momenta and frequencies: if a microscopic mechanism 
leads to a predominance of the forward scattering in 
the electron-phonon interaction, Tc is strongly enhanced. 
This effect was schematized in the past with the intro- 
duction of an effective cut-off in the electron-phonon in- 
teraction (Qc), which was argued to be due to strong 
correlations effects. 

In this work we have related the existence of Qc with 
the reduction of the screening properties due to correla- 
tion of a metal approaching a metal-insulator transition. 
The same effects which are responsible for the reduction 
of the screening (namely the loss of k-space coherence) 
are however also strongly detrimental to superconductiv- 
ity. In this work we have analyzed how the interplay 
between these effects is reflected on a Tc vs doping phase 
diagram. 

We have introduced a simple analytical model to sim- 
ulate the effects of the strong correlation on the one elec- 
tron Green's function. This model has also been em- 
ployed to estimate the role of the electronic screening 
on the electron-phonon scattering in correlated systems. 
We have shown that the reduction of the metallic charac- 
ter due to the electronic correlation implies a reduction 
of the "effective" Thomas-Fermi screening approaching 
6 — 0, where correlation is stronger. This results in a 
predominance of forward (small-q) scattering, which has 
been parametrized in terms of a phonon momentum cut- 
off Qc = fcTF/2A;F, where fcxr and are respectively the 
Thomas-Fermi and the Fermi vectors. 

We have also shown how the different parts of the 
electronic Green's function contribute to the supercon- 
ducting pairing; in particular, we have shown that the 
superconducting critical temperature is mainly deter- 
mined by the coherent excitations. The similarities 
and differences between our model and the two-band 
superconductivit y^^i^^ have also been discussed. 

Solving the nonadiabatic generalized ME equations, we 
obtained a Tc vs. 6 diagram, which can be ideally divided 
into three regions: 

(a) a high doping region, where superconductivity is 
suppressed by the negative contribution of the 
nonadiabatic channels to the electron-phonon pair- 
ing; 

(6) an extremely low doping region, where the poor 
metallic character is reflected in a vanishing coher- 
ent spectral weight. In this region, superconductiv- 
ity is extremely unstable and it can be overwhelmed 
by other electronic or structural instabilities in- 
duced by spin and charge degrees of freedom (anti- 
ferromagnetic fluctuations, stripes, chargc-density- 
waves, pseudogaps, . . . ). 
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(c) an intermediate doping region, in which the loss of 
coherent spectral weight is not large enough to pre- 
vent superconductivity, which is in turn enhanced 
by the positive contribution of the nonadiabatic 
channels of interactions. 

The resulting phase diagram bares strong resemblance 
with that of cuprates. We have in fact an overdoped re- 
gion, where superconductivity is triggered on by the pos- 
itive contribution of the nonadiabatic channels as doping 
is decreased; an optimal doping, where the enhancement 
due to the nonadiabatic interaction is counterbalanced 
by the reduction of the metallic character, and an under- 
doped region, where superconductivity disappears due to 
the incipient metal-insulator transition. In the qualita- 
tive scenario outlined here the origin of superconductivity 
in cuprates can be understood by focusing on the over- 
doped region, where the materials retain defined metallic 
properties; on the other hand, the exotic phenomenology 
of the underdoped region is only marginal. The occur- 
rence of different kinds of electronic/structural instabili- 
ties, not discussed in the present paper, is thus thought 
to be a by-product of the loss of metallic character which 
also drives the suppression of Tc as 5 more than to 
be the secret of the superconducting pairing. 

Once more, we wish to stress that what we present in 
this paper is a general scenario, based on the microscop- 
ical description of the interplay between nonadiabatic ef- 
fects and strong electronic correlation. A quantitative 
understanding of the specific phase diagram of cuprates 
should of course take into account specific features of 
these materials, such as Van Hove singularities and the 
d-wave symmetry of the gap. The possibility of a d-wave 
pairing within the context of electron-phonon supercon- 
ductivity has been discussed elsewhere^; we remind here 
that c?-wave symmetry was shown to be favoured by 
forward scattering, which in our model is enhanced as 
5^0. 

The qualitative behaviour of our results would be how- 
ever left unchanged by the inclusion of these effects. 

On the experimental ground we observe that the 
present scenario is supported by a detailed analysis of 
Tc vs. normal state properties in different families of 
cuprates. In Ref. 0, for example, the complex behaviour 
of Tc approaching the metal-insulator transition either 
by reducing the doping or by increasing the disorder was 
nicely pointed out by Osofsky et ai. The relation be- 
tween Tc and reduced screening properties was also dis- 
cussed there. Although we do not attempt to discuss 
the scaling relations close to the metal-insulator transi- 
tion in region (6), where a more specific treatment of the 
electronic correlation is needed, we think our analysis is 
somehow complementary to that of Ref. ^3- This sce- 
nario can also open new perspectives on the remarkable 
increase of Tc in granular metals and alloys4i 

Furthermore, a strong doping dependence of the 
electron-phonon properties in cuprates has also been re- 
ported by inelastic X-ray measurements of the phonon 
dispersion.® Experimental data in NCCO compounds 



were shown to be compatible with the theoretical calcu- 
lations, based on the shell model, assuming a negligible 
Thomas-Fermi vector {Qc = 0) for the strongest corre- 
lated undoped compound {S — 0), whereas a Thomas- 
Fermi momentum /ctf — 0.39A^^, comparable to that 
for LSCO, was estimated for 6 ~ 0.14. The corre- 
sponding dimensionless cut-off Qc would be hence es- 
timated Qc — 0.26 by using an in-plane Fermi vector 
kf ~ 0.741-1. 

L.B. wishes to thank Massimo Capone and Giorgio 
Sangiovanni for useful discussion and comments. The 
authors acknowledge financial support by the MIUR 
projects COFIN2001 and FIRB-RBAU017S8R, and by 
the INFM project PRA-UMBRA. 



APPENDIX A: GUTZWILLER SOLUTION FOR 
GENERIC U AND n 

In this appendix we provide a brief overview about the 
analytical solution of the Gutzwiller approximation for 
generic filling and Hubbard repulsion. 

Let us write the Hubbard Hamiltonian within the 
Gutzwiller approximation (in the paramagnetic state) as: 



H 



--fiU,n,d)\e\ + Ud, 



(AI) 



where n is the total electron filling, d the density of dou- 
ble occupancy sites, e the kinetic energy for site, and 



7(J7, n, d) = 



2{n-2d) 



\/l + + Vd 



n{2 — n) 

Minimizing Eq. ljAl|l with respect to d yields: 



= 



Vl - n 
U 



+2 



Uc 



•rd + Vd 
n{2 — n), 



(A2) 



Vd 



y^d{l -n + d) 

(A3) 



where we have introduced as usual Uc — 8|e|. 

After expanding the squares + ^TTT] ^ one can 

now isolate on the right side the remaining square roots: 



U 



2n{2-n) ( — 



(1 - 2n 4- Ad) 



and, by squaring both the sides of Eq. ljA4|l . all the 
remaining square roots are removed and we are left with 
a third order polynomial expression for d. We obtain 
namely: 



A^d^ + + Aid -F Ao = 0, 



(A5) 
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where 

^3 = 

A2 - 



16n(2-n)(^ 
M2-n)(^ 



(A6) 



U 



n{2-n)\ — ] -6n + 5(|A7) 



Ai = {l-n) 



+4n(2 - n)(l - 2n) [ — ] - n 



U 



An 



,(A8) 



(A9) 



Eq. I)A5|I can be easily solved to obtain dmin, and, in 
the standard notations, the Gutzwiller factor Z(U, n) = 
7(t/,n,d„ii„). 



APPENDIX B: ANALYTICAL EXPRESSION OF 
DIFFERENT PHYSICAL QUANTITIES 

1. Thomas- Fermi Screening 



In this section we provide some useful analytical ex- 
pressions for the different contributions (IIc-c, Hc-i, 
ni_i) to the response function 11 involved in the eval- 
uation of the Thomas-Fermi screening as limit fc^p cx 
limq_^on(q,w = 0). 

In the RPA approximation n(q, oj) is given by: 

n(q,cj) = -|-^ / duj'G{k + q,uj + Lj')Gik,uj'). (Bl) 

We employ the simple model of Eqs. l@J, ((SJ for the 
electron Green's function in the presence of correlation. 

From simple scaling relation it is straightforward to 
recognize that 



ZIl{Z = 0) 



-ZNo. 



(B2) 



The analytical expressions for Hc-i, Ili^i are straight- 
forward but more cumbersome since they involved the 
explicit integration over the upper/lower Hubbard bands. 
One obtains: 



{u/2 



^-{l-^)E/2) X 



X In (^C//2 + - (1 - ■^)E/2 

+ {U/2 + nE/A + ZE/2) x 
X \n{U/2 + nE/4: + ZE/2) + 

- {U/2-nE/A + ZE/2) x 

X \TL{U/2-nE/A + ZE/2) + 

- {U/2 + nE/A - ^i) X 

X \TL{U/2-nE/A- ii) + 
+ {U/2 - nE/A + ti)x 

X \n{U/2-nE/A + ^i) (B3) 

The incoherent-incoherent contribution gives: 

Bi-i = -2N^{l-Zf[{U + E/2)ln{U + E/2)+ 
+ {U - E/2)\n{U ~ E/2) + 

- [U +{n-l)E/2)\n{U +{n-l)E/2) + 

- [U +{l-n)E/2)\n{U +{l-n)E/2)] 

2. Superconducting properties 

Here we report the explicit expressions for the coherent 
contribution to different quantities in Eqs. I|19|l - (|2()(l . 

Let us consider for instance rjm[Gcoh\- In this case 
?7m[Gcoh] = VmlGcoh] and: 



'7m[Gcoh] = / deN{€) 
Z 



Euj„ 



iuj„i — Zei^ — /i —iuj. 
fZE - fi 



arctan 



arctan 



fZE+ a 



The expression (IB4|I corresponds just to the i]m{E) for 
an uncorrelated system with reduced spectral weight Z 
and rescaled bandwidth ZE: ri,n[Gcoh]{E) = Zrim{ZE). 
Similar considerations apply for the vertex and cross 
function: P[Gcoh](-E', Qc; n^) — ZP{ZE,Qc]n,ra), 
C[Gcoh\{.E,Q(.;n,m) ~ ZC{ZE,Qc;n,m), where 
P{E, Qc, n, m) and P{E, Qc] n, m) in the absence of elec- 
tronic correlation were computed in Refs. .17Ml8ii36i : 

P{E, g,; n, m) = D(w„ - w/)|B(r^, m, + 



H, 



2N{{)f{\-Z) (^ZE/2 + U/2 + {l-^)E/2^ X ^ A{n, to, I) ~ B{n, m, - uji-n+m) 



X \ni^ZE/2 + U/2+ {1- -)E/2 



ZE/2 + U/2 - (1 - ■^)E/2j X 
X \tl{zE/2 + U/2- {l-^)E/2}j 
^ (c//2 + M+(l-|)£;/2) X 
X ln([//2 + ^+(l-|)i;/2) 




(B5) 



C{E, Qc] n, m) = -D(w„ - uJi)D{uoi - curn) x 
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X < 2B{n, — m, I) + arctan 



2EQI 



— arctan 



\^i — <^l-n+m\ 
A{n, -m, I) - B{n, -m, l){uji - oji-n-m)^ 

EQl\u)l - Wl-n-m\ 

where 

A{n,m,l) = {UJI - UJl-n+m) ^^^^^"^ + 



E 



>iB6) 



B{n,m,l) = {uJi~Ui-n+m) 



EcJi- 



2[{E/2r+wl 



(B7) 

+ 



■n+mj 



E 



2[{E/2)^ + ul^^J 



(B8) 



^ M.K. Crawford, W.E. Farneth, E.M. McCarron, R.L. Har- 
low, and E.H. Moudden, Science 250, 1390 (1990). 

^ J.P. Franck, S. Gygax, S. Socrcnscn, E. Altshuler, A. 
Hnatiw, J. Jang, M.A.-K. Mohamed, M.K. Yu, G.I. 
Sproule, J. Chrzanowski, and J.C Irwin, Physica C 185- 
189, 1379 (1991). 

^ See H. Keller, Physica 326, 283 (2003) and references 
therein. 

" F. Raffa, T. Ohno, M. Mali, J. Roos, D. Brinkmann, K. 

Conder, and M. Eremin, Phys. Rev. Lett. 81, 5912 (1998). 
^ D. Rubio Temprano, J. Mesot, S. Janssen, K. Conder, A. 

Furrer, H. Mutka, and K. A. Mller, Phys. Rev. Lett. 84, 

1990 (2000). 

® M. d'Astuto, P. K. Mang, P. Giura, A. Shukla, P. Ghigna, 

A. Mironc, M. Bradcn, M. Graven, M. Krisch, and F. Sette, 
Phys. Rev. Lett. 88, 167002 (2002). 

B. Friedl, C. Thomsen and M. Cardona, Phys. Rev. Lett. 
65, 915 (1990). 

* A. Lanzara, P.V. Bogdanov, X.J. Zhou, S.A. Kellar, D.L. 

Feng, E.D. Lu, T. Yoshida, H. Eisaki, A. Fujimori, J.-I. 

Shimoyama, T. Noda, S. Uchida, Z. Hussain, and Z.-X. 

Shen, Nature 412, 510 (2001). 
^ M. L. Kulic, Phys. Rep. 338, 1 (2000). 
1° R. Zeyher and M.L. Kulic, Phys. Rev. B 53, 2850 (1996). 
" M. Grilh and C. Castellani, Phys. Rev. B 50, 16880 (1994). 
^'^ Z.B. Huang, W. Hanke, E. Arrigoni, and D.J. Scalapino, 

cond-mat/030613 (2003). 

M. Wcgcr, M. Potor, and L.P.Pitacvskii, Z. Phys. B 101, 
573 (1996). 

" D. Fay and M. Weger, Phys. Rev. B 62, 15208 (2000). 

A.A. Abrikosov, Physica C 222, 191 (1994). 
^'^ R. Chitra and G. Kothar, Phys. Rev. Lett. 84, 3678, 

(2000). 

^''^ L. Pietronero, S. Strassler and C. Grimaldi, Phys. Rev. B 

52, 10516 (1995) 
^® C. Grimaldi, L. Pietronero, S. Strassler, Phys. Rev. B 52, 

10530 (1995). 

C. Grimaldi, L. Pietronero, S. Strassler, Phys. Rev. Lett. 

75, 1158 (1995). 
^° A. B. Migdal, Zh. Eksp. Teor. Fiz. 34, 1438 (1958) [Sov. 



Phys. JETP 7, 996 (1958)]. 

A.I. Liechtenstein and M.L. Kulic, Physica C 245, 186 
(1995). 

22 A.A. Abrikosov, Physica C 214, 107 (1993). 

23 J. Bouvier and J. Bok, Physica C 249, 117 (1995). 

2'' I. Chang, J. Friedel, and M. Kohmoto, Europhys. Lett. 50, 
782 (2000). 

P. Paci, C. Grimaldi, L. Pietronero, Eur. Phys. Journ. B 
17, 235 (2000). 

2^ A.A. Abrikosov, L.P. Gorkov, and I.E. Dzyaloshinski, 

Methods of Quantum Field Theory in Statistical Phisics 

(Dover, New York, 1963). 
2^^ For a review: Correlated Electron Systems, ed. by V.J. 

Emery (World Scientific, Singapore, 1993). 
^® J. Hubbard, Proc. R. Soc. London A 276, 238 (1963). 
2® M.C. Gutzwiller, Phys. Rev. Lett. 10, 159 (1963). 
3° G. Kotliar and A.E. Ruckenstein, Phys. Rev. Lett. 57, 1362 

(1986). 

A. Georges, G. Kotliar, W. Krauth, and M.J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

P.B. Allen and B. Mitrovic, in Solid State Physics, v. 37. 
ed. by H. Ehrcnrcich, D. TurnbuU, and F. Seitz (Academic, 
New York, 1982) 

33 G.D. Mahan, Many-Particle Physics 3rd Edition (Kluwer 
Academic, New York, 2000). 

34 W.E. Pickett, Rev. Mod. Phys. 61, 43 (1989). 

35 G. Baym and L. Kadanoff, Phys. Rev. 124, 287 (1961). 
3^ M. Scattoni, C. Grimaldi and L. Pietronero, Europhys. 

Lett. 47, 588 (1999). 
3^" E. Cappelluti, C. Grimaldi and L. Pietronero, Phys. Rev. 

B 64, 125104 (2001). 
3® H. Suhl et al, Phys. Rev. Lett. 3, 552 (1959). 
3^ A.Y. Liu, I.I. Mazin and J. Kortus, Phys. Rev. Lett. 87, 

087005 (2001). 

4° M.S. Osofsky, R.J. Soulcn, Jr., J.H. Classen, G. Trotter, 
H. Kim, and J. Horwitz, Phys. Rev. B 66, 020502 (2002). 
M.S. Osofsky, R.J. Soulen, Jr., J.H. Classen, G. Trotter, H. 
Kim, and J. Horwitz, Phys. Rev. Lett. 87, 197004 (2001). 



